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We provide an analytic solution to the short wave length limit of the integro-differential equation 
describing the damping of the tensor modes of gravitational waves. 
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I/"") ' In a recent paper Q], Weinberg gives an integro-differential equation for the propagation of cosmological 

gravitational waves. In particular he writes an equation for the perturbation to the metric hij(x.,t) and then 
' defines x( u ) as 

(o)x(«), (i) 

\ where u is the conformal time multiplied by the wavenumber 

, f* df 

u = k W)- (2) 

x( u ) satisfies an integro-differential equation which for short wavelengths (wavelengths which entered the horizon 
■ while the universe was still radiation dominated) is given by P, Q 

O ■ r 
0\ u 2 X "(u) + 2u X '(u) + u 2 X {u) = -24/„(0) / dU K(u - U) X '(U) . (3) 

O Jo 

in 

The fraction of the energy density in neutrinos is f v (0) = 0.40523 and the kernel K will be discussed in detail 
below. The initial conditions are 



6 



X (0) = 1, X'(0) = 0. (4) 



In the absence of free-streaming neutrinos the right-hand-side of Eq.JSJ) is zero and x( u ) = sin(u)/ti. The 
suppression of these modes is due to the presence of the neutrinos where the solution of Eq. (0 approaches, for 
u » 1, 



x{u) — ► A sin(u + S)/u , (5) 



and the value of A 2 is the quantitative measure of that suppression. After deriving the above Weinberg says 
"A numerical solution of Eqs. (22) and (23) [our Eqs.© and shows that x( u ) follows the /„ = solution 
pretty accurately until u ks 1, when the perturbation enters the horizon, and thereafter rapidly approaches the 
asymptotic form (Eq.JSJ)), with A = 0.8026 and 5 very small." 

Our purpose here is to provide an analytic solution of Eqs.|(2J) and Jlj so that readers can verify the quoted 
statement for themselves. 

The importance of these results is shown by the remainder of the quote "This asymptotic form provides the 
initial condition for the later period when the matter energy density becomes first comparable to and then 
greater than that of radiation, so the effect of neutrino damping at these later times is still only to reduce the 
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tensor amplitude by the same factor A = 0.8026. Hence, for wavelengths that enter the horizon after electron- 
positron annihilation and well before radiation-matter equality, all quadratic effects of the tensor modes in the 
cosmic microwave background, such as the tensor contribution to the temperature multipole coefficients Ce and 
the whole of the "B-B" polarization multipole coefficients Ces, are 35.6% less than they would be without the 
damping due to free-streaming neutrinos." 

A solution to Eq. J2Jl is a series of spherical Bessel functions Q 

oo 

x{u) = ^2 a n j n (u) . (6) 

Inserting Eq. JBJ in the left-hand-side of Eq. and using the differential equation for spherical Bessel functions 
leaves 

oo 

n(n + l)a n j n (u) . (7) 

n=0 

The right-hand-side of Eq. J^J requires more work. The kernel is itself a sum of spherical Bessel functions 

K(u) = 1 J dx(l-x 2 ) 2 e" u (8) 

sin u 3 cos u 3 sin u , , 

= 5 Z— + 9 ) 

U 3 U 4 M 5 W 

~\- f , . 10 . . , 3 . A 

Jo(u) + —J2(u) + -u(u) , (10) 



15 V 27" v ' V 

and the derivative of x( u ) is given by 

oo 

X '(U) = J2 a nfn(u) (11) 

n=0 

V- N„_i(m) - (n + l)j n+1 {u)} 
So the RHS of Eq. © is -CI(u) with C = 1.6 /„(0) = 0.648368 and I(u) given by 

OO 

= d m ^ " I n , m {u) (13) 

m=o,2,4 n=o ^ n + ^ 

where 

/n,m(«) = / dl/ Jm(« ~ U) [n Jn -l{U) - (fl + l)j„+l(U)] . (14) 



The d m are given by Eq. (|1()|> where we have factored out a j^. Evaluating I n ,m is an exercise in using Abramowitz 
and Stegun [J|(AS). First, use the fact that the Fourier transform of a Legendre polynomial is a spherical Bessel 
function (AS 10.1.14) 

(-i) n r 1 

Mx)= [ —^ dse lxs P n {s) (15) 



" J -1 

to replace both Bessel functions in Eq. I|14(l . This makes the integral over U trivial and we have 

In.rn(u) = ^ / ds / dt P m (s)[nP n ^(t) + (n + l)P„+i(t)] . (16) 

4 J -i J -i t-s 



Now use the definition of the Lcgcndrc function of the second kind, 

1 f 1 

Qn(z) = -J dxiz-X^P^x), (17) 

(AS 8.8.3), to evaluate the integral in Eq. (|16[) over the variable that does not appear in the exponent (s in the 
first term, t in the second) to obtain 

I n ,m(u) = ^-L J dte ltu {Q m (t)[nP n ^(t) + {n + l)P„+i(t)] + P ro (f)[nQ„_i(t) + (n + l)Q n+1 (t)]} . (18) 

Next, by replacing the remaining exponential with the familiar expression from quantum mechanical scattering, 

e uu = Y^(2i+l)i e Mu)Pi(t), (19) 

£ 

the expression for I n ,m{u) becomes 

I n , m {u) = J2 ^r^ HT +m ~ e M(u>) { f dtP e (t)Q m (t)[nP n ^(t) + (n + l)P n+1 (t)] 



l 

r-l 

dtP t (t)P m (t)[nQ n -i(t) + (n + l)Q„ +1 (t)] \ . 



Eq. lO can be simplified using (AS 8.6.19), 



where 



and the formula 



,(a:) = \p m {x) In ~ W m _i(x) , (21) 

2 1 — x 



, . ^ v— ^ ' 2m — 4fc — 1 _ 
W m ^(x)= g ^I^WiW, (22) 



P/(x)P m (x) = |<^,0,m,0|L,0>| 2 P L (x) (23) 

L=\t-rn\ 



to express p?(x)P m (x) in terms of Pj,(x)'s and, with the aid of Eq. I|21|l . Pi{x)Q m {x) in terms of <5l(x)'s as 



Pi(x)Q m (x)= J2 [\<t,0,rn,0\L,0>\ 2 (Q L (x)+W L - 1 (x))]~P e (x)W m - 1 (x). (24) 

L=\l-m\ 

Finally, the terms in Eq. ij2"U|) involving the products of P n (x)'s and Q m (x)'s cancel using (AS 8.14.10), 

J dx(Q L (x)P m±1 (x)+P L (x)Q m±1 (x))=0, (25) 

and I n . m (u) reduces to 

2/4-1 r f 1 £+m 

In, m {u) = H)" +TO "^N{ / dt J2 I < ^, 0, m, 0|i,0 > \ 2 WL-i(i) [nP n -i(t) + (n+ l)P n+ i(t)] 

L=\i-m\ 

^ dtP e (t)W m ^(t)[nP n ^(t) + (n + l)P„ +1 (t)]} (26) 
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The contributions to the coefficient of each ji{u) in Eq. i|26|) can be straight-forwardly evaluated; those in the 
sum by directly using orthogonality and the remaining terms by expressing the product of two P/s as a sum 
of Pi 's and again using orthogonality. The orthogonality of the Legendre functions means that the I which is 
summed over in Eq. I|26|l can only take on the values n + 2k where k = 0, 1, 2, • • • so we replace 

m=0,2A 



with 



k=0 



i.e., the sum over I in I n ,m{u) is replaced by a sum over k and each c n ^ is the sum of the contributions from the 
three terms in the kernel, m = 0, 2, 4. Setting Eq. Q equal to — CI(u) we have 

oo oo 

^ n(n + l)a n j n (u) = -C a n c n ,kjn+2k{u) , (27) 

n=0 n,k=0 

where the c n ,k are known numbers and we can find the expansion coefficients, a n , recursively by equating the 
coefficients of each order Bessel function in Eq. (|27(l . 

The coefficients of ji (it) in Eq. l|27|l give 

2cn = -Cci.001, (28) 

where c„.o is equal to 1 — 5 n ,o so c^o is 1. The only solution of this equation is a\ = 0. This ensures that 
the second of the initial conditions, Eq. Q, is satisfied. The equality of the coefficient of j'3 (it) shows that 03 
is proportional to a±. Similarly a§ is a linear combination of a\ and 03, 07 a linear combination of ai, 03, and 
05, • • •. Thus the coefficients of all the odd order Bessel functions in Eq.© are zero. There is no mixing between 
the coefficients of the odd order Bessel functions and those of even order because the Clcbsch-Gordan coefficients 
< a, 0, b, 0|c, > are zero if a + b + c is an odd number. 
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TABLE I: The coefficients needed for Eq. IpTJ to evaluate Eq. © up to n = 12. 



Thus the only nonzero a n in Eq. JSJ are those with even n. a Q does not appear in Eq. I|27|) but is determined 
by the first of the initial conditions of Eq. which fixes it to be unity. The c Ui k necessary to find 02, • • • , «i2 
are shown in Table [I] 
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The equations for these a n can be read off from Eq. 1)2 7)1. Using ao = 1, we have 

«2 = — O ■ 



0,4 = — C 



6 + C 

CQ,4 + Q2C2,2 

20 + C 

(29) 



G2ri 



n-1 

C ^ a2feC2fe,2n-2fc 
fc=0 

2n(2n+l) + C 
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FIG. 1: (color online) The solution x( u ) to Eq. © (solid) is compared to the f v = solution ji'o(w) (dotdashed) and the 
asymptotic solution 0.80313 jo (it) (dashed). 

Numerically wc find that the a n decrease very quickly 



"0 


= 1.0 


a 2 


= 0.243807 


C<4 


= 5.28424 x 10" 2 


066 


= 6.13545 x 10~ 3 


<"<8 


= 2.97534 x 10~ 4 


aio 


= 6.16273 x 10" 5 


ai2 


= -4.78885 x 10" 6 



For large argument, all of the even order Bessel functions go as ± sin x/x so the A in Eq. JSJ) is 

5 

A = Et- 1 )" 02 ™' 

71 = 

= 0.80313. (30) 
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Since there are no odd order Bessel functions in the expansion, the phase S in Eq. JSJ) is zero. 

The extent to which xi u ) departs from the the /„(0) = solution jo(u) and approaches the asymptotic solution 
0.80313 jo(u) is illustrated in Fig^We see that x( u ) is discernably different from jo(u) by u = 1.5 — 2 and that 
it is essentially identical to 0.80313 jo(u) by u = 5. 
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